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ABSTRACT 

Since 2002, the number of detected blazars at gamma-ray energies above 100 GeV has more 
than doubled. I study 17 blazars currently known to emit E > 100 GeV gamma rays. Their 
intrinsic energy spectra are reconstructed by removing extragalactic background light atten- 
uation effects. Luminosity and spectral slope in the E > 100 GeV region are then compared 
and correlated among each other, with X-ray, optical and radio data, and with the estimated 
black hole (BH) masses of the respective host galaxies. 

According to expectations from synchrotron self-Compton emission models, a correlation 
on the 3.6-cr significance level between gamma-ray and X-ray fluxes is found, while corre- 
lations between gamma-ray and optical/radio fluxes are less pronounced. Further, a general 
hardening of the blazar spectra in the £ > 100 GeV region with increasing gamma-ray lumi- 
nosity is observed, both for the full 17-source sample and for those sources which have been 
detected at distinct flux levels. This goes in line with a correlation of the gamma-ray luminos- 
ity and the synchrotron peak frequency, which is also seen. Tests for possible selection effects 
reveal a hardening of the spectra with increasing redshift. The blazar gamma-ray emission 
might depend on the mass of the central BH. The blazars under study do, however, show no 
correlation of the BH masses with the spectral index and the luminosity in the £ > 100 GeV 
region. 

I also consider temporal properties of the X-ray and £ > 100 GeV gamma-ray flux. No 
general trends are found, except for the observation that the blazars with the most massive BHs 
do not show particularly high duty cycles. These blazars include MknSOl and PKS 2155-304, 
for which recently very fast flares have been reported. In general, VHE flare time-scales are 
not found to scale with the BH mass. 

As a specific application of the luminosity study, a constraint for the still undetermined 
redshift of the blazar PG 1553+113 is discussed. 

Key words: galaxies: active - BL Lacertae objects: individual (lES 0229-1-200, lES 0347- 
121, IBS 1011+496, lES 1101-232, lES 1218+304, IBS 1959+650, lES 2344+514, 3C279, 
BL Lacertae, H 1426+428, H2356-309, Mknl80, Mkn421, Mkn501, PG 1553+113, 
PKS 0548-322, PKS 2005-489, PKS 2155-304) - black hole physics - galaxies: jets. 



1 INTRODUCTION 

All but one of the detected extragalactic very high energy (VHE, 
defined by £ > 100 GeV) gamm a (y) ray sources so far a re blazars. 
Within the unified scheme (e.g. lUrrv & Padovanilll995h of active 
galactic nuclei (AGN), blazars comprise the rare and extreme sub- 
classes of BL Lac objects and flat spectrum radio quasars (FSRQs). 
These are characterised by high apparent luminosities, short vari- 
ability time-scales, and apparent superluminal motion of jet com- 
ponents. These observations can be explained by highly relativis- 
tic, beamed plasma outflows (jets) close ly aligned to the observer's 
line of sight JBlandford & Koni gill 19791) powered by central super- 
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massive black hol es accreting at sub-Eddington rates iLvnden-Belll 
[l969l; |Rees||l97^b). The prime scientific interest in VHE y-ray 
emitting blazars (in the following, 'VHE blazars') is twofold: (1) 
To understand the particle acceleration and y-ray production mech- 
anisms, assumed to take place in the jets and to be linked to the cen- 
tral supermassive black hole (BH). Knowledge of the VHE emis- 
sion process will also contribute to the further understanding of the 
accretion processes in AGN, jet formation processes, and the jet 
structure. (2) To use the VHE y-rays as a p robe of the extragalac- 
tic background light (EBL; e.g. lHauser & D wek 2001; KashlinskyI 
1 20051) spectrum in the wavelength range between about 0.3 to 
30 fim. Determining the EBL spectrum in this wavelength range 
may allow to constrain the star formation rate (convolved with the 
initial mass function) in the early Universe. In order to assess both 
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Figure 1. Currently known VHE y-ray blazars along with the identified (66 objects) and tentatively identified AGN (27 objects) in the 3rd EGRET catalogue 
of y-ray sources (solid grey dots: identified AGN; open grey dots: tentatively identified AGN). EGRET data from .Hartman et aL 1 ,1999. ). The sources are shown 
in a galactic coordinate system. 



issues, it is essential to have a large sample of VHE y-ray blazars at 
hand. Ideally it should encompass a wide range in redshift for EBL 
studies and at the same time include groups of sources at similar 
distances in order to probe and compare properties of the individ- 
ual sources without possible systematic uncertainties caused by the 
EBL de-absorption. 

The preconditions for such studies have much improved re- 
cently: Before 2004, only a few nearby extr agalactic so urces had 
been established as VHE y-ray emitters (e.g. Morill2003l) and pro- 
vided hardly enough data to perform comparative studies. Around 
2004, the third gene ration of imaging air Cerenkov telesco pes 
(lACTs, e.g. MAGIC. iBaixeras et al.ll2004ICortina et al.ll2005l and 



H.E.S.S.. lHintonll2004l) . the most successful tools so far to explore 
VHE y-rays, started to deliver scientific results. To date, the VHE 
blazar sample with available spectral information in the VHE re- 
gion comprises 17 BL Lac objects, among them one LBL object, 
BL Lacertae. Furthermore, now the redshifts of the known VHE 
blazars reach up to z = 0.212 - or even to z = 0.536, considering 
the recently announced discovery of the first FSRQ in VHE y-rays, 
3C 279 (Teshima et al. 2007). For 3C 279, however, no VHE spec- 
trum has been published yet, therefore the VHE luminosity and 
spectral slope of 3C 279 are not yet available for this study. M87, 
a FR I radio galaxy also detected in VHE y-rays, is not included in 
the study, as its VHE y-ray prod uction mechanism may difl"er from 
that in blazars JAharonian et al.i 2006c ). Fig.[T]shows the sky posi- 
tions of the known VHE blazars in galactic coordinates along with 
the AGN identified in the 3rd EGRET catalogue. 

The electromagnetic continuum spectra of blazars extend over 
many orders of magnitude from radio frequencies to sometimes 
multi-TeV energies and are dominated by non-thermal emission 
that consists in a vF,, representation of two pronounced peaks. 
The low-energy peak, located between the IR and hard X-rays, is 
thought to arise from synchrotron emission of ultrarelativistic elec- 
trons, accelerated by shocks moving along the jets at relativistic 
bulk speed. Depending on the location of the low-energy peak, 
BL Lac objects are often referred to as high-frequency peaked 
(HBL; in the UV to X-ray domain) or low- frequency peaked (L BL; 
in the near-IR to optical) BL Lac objects dFossati et al.ll 19981) . al- 
though the transition is smooth rather than dichotomic. The ori- 
gin of the high-energy peak at MeV to TeV energies is still de- 
bated. It is commonly explained by inverse Compton upscatter- 
ing of low-energy photons by electrons. The seed photons may 



originate from synchrotron radiation produced by the same elec- 
tron population (synchrotron-self Compton (SSC ) models; e.g. 
iMaraschi, Ghisellini & Celottil 1 19921 ; ICoppil Il992h or belong to 
ambient thermal photon fields (external inverse Compton mod- 
els; e.g. iM elia & Konigl 1989; ISikora, Begelman & Reed Il994l ; 
IPermer & Schlickeiseriil994) . In hadronic models, which can also 
explain the observed features, interactions of a highly relativistic je t 
outflow with ambient matter (Dar & Laor ' 19971 ; iBednarekll 1993h . 
proton-indu ced cascades (Mannheim 1993). svnchrotro n radiation 
by protons ( lAharoniarj|200d ; lMucke & Protherodl 20011) . or curva- 
ture radiation, are responsible for the high energy photons. The 
AGN identified in the EGRET data are predominantly powerful FS- 
RQs and quasars with SEDs peaking at rather low frequencies, and 
thus only few of these (Mkn 421, PKS 2155-304, BL Lacertae and 
3C 279) were also detected in the VHE range. 

Knowing the variability time-scales of the VHE y-ray emis- 
sion and the form of the two-bump spectral energy distribution 
(SED) enable s the derivation of all input parameters of one-zone 
SSC models dXavecchio et alj|l99 8). which describe the observed 
emission in BL Lac objects reasonably well. Strictly simultaneous 
and temporally-resolved measurements of the SED, however, are 
only rarely possible and more often than not are also severely re- 
stricted by the (temporal) instrumental resolution. Generally, de- 
tailed spectral studies particularly in low-emission states are rather 
demanding. In addition, the determination of the location of the 
high-energy peak v'^^ generally would require complementary 
satellite detector coverage of the SED between some hundred MeV 
and X 50 GeV This region of the SED, however, is difficult to ac- 
cess due to the low fluxes expected (up to » 100 MeV) from extra- 
galactic sources in between the two bumps and due to insufficient 
instrumental sensitivity for energies exceeding some GeVs. 

In this paper, for the first time studies of the VHE emission 
properties of the complete set of all currently known VHE blazars 
(cf. Tab. [T} are performed. First, the detected VHE blazars are 
brought into context with the AGN searches conducted by lACTs 
so far and the expected y-ray attenuation by the EBL in Sect. [2] Af- 
ter a study of the black hole mass distribution of the VHE blazars 
in Sect. [3] I infer intrinsic emission properties in the VHE y-ray 
regime in Sect.|4] Because measurements of the location and shape 
of the high-energy bump are elusive at present for almost all VHE 
blazars, the observed y-ray luminosity and spectral slope in the 
VHE region are used as auxiliary observables to characterise the 
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Table 1. Extragalactic VHE y-ray sources, listed in chronological order of 
their discovery. 



Source 


Type 


Redshift z 


Discovery reference 


Mkn421 


HBL 


0.030 


Punch et al. (1992) 


MknSOl 


HBL 


0.034 


Ouinn et al. (1996) 


lES 2344+514 


HBL 


0.044 


Catanese et al. (1998) 


lES 1959+650 


HBL 


0.047 


Nishivamaetal. (1999) 


PKS 2155-304 


HBL 


0.116 


Chadwick et al._ (1999) 


H 1426+428 


HBL 


0.129 


Horan et al. (2002) 


M87° 


FRl 


0.0044 


Aharonian et al. (2004) 






PKS 2005-489 


HBL 
HBL 


0.071 
0.182 


Aharonian et al. (2005b) 


lES 1218+304 


Albert et al. (2006b) 


H 2356-309 


HBL 


0.165 


Aharonian et al. (2006a) 


lES 1101-232 


HBL 


0.186 


Aharonian et al. (2006a) 


PG 1553+113 


HBL 


b 


Aharonian et al. (2006b). 




_Albert et al. (20()7a) 


Mknl80 


HBL 


0.045 


Albert et al. (20()6c) 


PKS 0548-322 


HBL 


0.069 


Superina et al. (2007) 


BL Lacertae 


LBL 


0.069 


Albert et al. (2007d) 


lES 1011+496 


HBL 


0.212 


Albert et al. (2007e) 


lES 0229+200 


HBL 


0.139 


Raue et al. (2007) 


lES 0347-121 


HBL 


0.188 


Aharonian et al. (2007c) 


3C 279" 


FSRQ 


0.536 


Teshima et al. (2007) 







The upper part of the table shows the confirmed sources prior to the advent 
of new generation instruments like MAGIC and H.E.S.S., while the lower 
panel summarises the sources discovered after 2002. HBL: High-frequency 
peaked BL Lac object, LBL: Low-frequency peaked BL Lac object, FSRQ: 
Flat spectrum radio quasar, FR: Fanaroft-Riley galaxy. " M87 is not in- 
cluded in the present study. This redshift is currently under discussion, cf. 
Sect. l5.6l " No spectrum of 3C 279 has been published yet. 



VHE y-ray emission. The main part of the paper (Sect. l5.1H5.5t 
is devoted to the search for correlations of these observables with 
the X-ray emission properties, the optical and the radio luminos- 
ity, and with the black hole mass estimations. In Sect. l5.6l the VHE 
blazar luminosity distribution is used to address the specific prob- 
lem of the unknown redshift of PG 1553+ 1 13 by deriving an upper 
redshift for this VHE blazar. Finally, Sect. 15.71 turns to the study 
of X-ray and VHE y-ray timing properties. Sect.[6]summarises the 
main conclusions of the studies. 



2 POPULATION STUDIES AND THE y-RAY HORIZON 

When travelling cosmological distances, VHE y-r ays interact with 



the low-energy photons of the EBL (see, e.g. iNi kishov 1962 



iGould & SchredeJl966l : [Hauser & Dw"e3l200ll : lKashlinskv..200S) . 
The predominant reaction y vhe + 7ebl — > e^ e" modifies source- 
intrinsic y-ray energy spectra. The cross-section of this process 
peaks strongly at iicM = 1.8x2meC", therefore a given VHE photon 
energy probes a narrow range of the EBL spectrum. The part of the 
EBL to which VHE y-rays are sensitive comprises the (redshifted) 
relic emission of galaxies and star-forming systems and the light 
absorbed and re-emitted by dust. The EBL attenuation results in a 
maximum distance over which photons wit h a particular energy can 
survive: The Fazio-Stecker relatio n (FSR, Ipazio & Steckeilll97d : 
IStecker. de Jager & Salamonlll992h describes the distance at which 
the optical depth for a VHE photon of a given energy reaches unity 
(attenuation by a factor e"' ). Thus the FSR defines the cosmological 
y-ray horizon. Fig. [2^ shows the instrumental low-energy thresh- 
olds from searches for VHE y-ray emission from AGN. Along with 



these, the FSR for different EBL models as given bv lKneiske et al.l 
(2004) is plotted. The models differ in the IR density, dust prop- 
erties and the star formation rate in the early Universe. The FSR 
divides the plot into a region from which no y-rays can reach the 
Earth and into another region, in which positive detections are to be 
expected or a too weak source-intrinsic emission made detections 
fail (due to insufficient instrumental sensitivity). Obviously, with 
a decreasing instrumental energy threshold, the visible Universe 
'opens up', providing access to a larger source population. CuiTent 
EBL models result in a steepening of the intrinsic spectra from « 
200 GeV on (power-law spectra are softened, but their shape is ap- 
proximately retained), while for lower energies the effects are min- 
imal. When an optical depth of one is reached, a quasi-exponential 
cutoff in the observed spectra occurs. Fig. [2j3 shows the energy 
ranges over which blazars in VHE y radiation have been detected. 
Up to now, only for the strong, close-by blazars Mkn 421 and 
Mkn 501 indications of the expected exponential high energy cut- 
off have been observed thanks to high y statistics JAharonian et al.l 
l2001al : lKrennrichetal.ll200lUAIbertetal.| l2007cV The observed 
spectra of all other blazars can be accurately described by power- 
laws or broken power-laws. While most of the nearby VHE blazars 
caimot constrain the current EBL models, some of the sources at 
z > 0. 1 start chall enging these jAharonian et alj|2006al . 12007 d : 
iTf.shima et al .1120070 . as no cutoff's have been observed so far at 
the high-energy ends of their y-ray spectra. 



3 BLACK HOLE MASSES IN BLAZARS 

It is well established that all galaxies with a massi ve bulge compo- 
nent host supermassive black hole s in their centres jRichstone et al.l 
119981 : iBender & Kormendv 120030 . There are a couple of indirect 
methods to infer the masses of the central BHs: One is to estimate 
M. using the correlation between M. and the central velocity dis- 
persion cr of the host galaxy (M, - cr relation. iFerrarese & MerrittI 
1200(1 : iGebhardt et al.l 12000) found from stellar and gas kinemat- 
ics and maser emission. I estimated the black hole masses of 
VHE y-ray em i tting b lazars using the M, - a relation given by 
iTremaine et alj ( l2002h . This approach assumes that AGN host 
galaxies are similar to non-active galaxies. The velocity dispersions 
were collecte d from the literature or ar e inferred from the funda- 
mental plane JDiorgovski & Davis! 19871) . a relation between cr, the 
effective galaxy radius R^, and the corresponding surface bright- 
ness (/Je), which is valid for elliptical galaxies, in particular also 
for AGN and radio galaxies JBettoni et al.ll200ll) . Whenever more 
than one cr value is given in the literature, individual masses were 
derived for each of the a values and averaged. The cr values and 
the resulting black hole masses for the VHE y-ray blazars studied 
here are given in Tab.|2] 

The determination of M. suffers from rather large sys- 
tematic uncertainties due to the different methods used to de- 
rive cr. The relation betwee n M. and bulge luminosity Lb 
dKormendv & Richstonelll995n has generally a larger scatter than 
the M. - cr relation and was therefore only used for PKS 2155- 
304, because for this blazar no cr or (yUe> measurement is available. 
I used the i?-band luminosity given by Falomo ( 1996) to calcu - 
late log(M./M(r) ) = 8.9 8 ± 0.44 using eq. 12 in lOrahanj ( l2007t) . 
lAharonian et al.l ilOOlih give an estimate of M. - (1...2) x 
IO'Mq. For a comparative study using values inferred by differ- 
ent methods is not advisable due to possible different systematics. 
Therefore, the M, value used for PKS 2155-304 should be taken 
with care in the following. The BH mass of 3C 279 was determined 
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Figure 2. (a) Lower energy thiesholds from searches for VHE y-ray emission from AGN. Dots: blazars, crosses: other AGN typ es (starburst galaxies , radio 
galaxies, etc.). The curves represent Faz io-Stecker relations (flux attenuation by a factor e~*) for diff erent EBL models given bv^Kneis ke et al, 120041). Data 
from searches by HEGRAJAharonian et aljl2004l) . H.E.S.S. ^haronian et al. 2005c; Benbow & Bi ihler 2007), Whipple I Kerrick et alJ ll995l : lHoranet alj 
|2004) . VERITAS jCogaijIJOO?! : lKrawczvnskill2007t) . CANGAROO iNishijima 2002) and MAGIC I' Albert et aUl200ai) . (b) Detected VHE y-ray sources (see 
Tab.|3]for references). Shown are the energy ranges of the measured y-ray spectra; high-energy cutoft's were found only for Mkn 421 andMknSOl, for all other 
sources arrows represent spectra possibly continuing to higher energies. The dashed lines represent PG 1553+113 at the lower and upper limit for its redshift, 
Z > 0.09 >Sbarufatti et aU20Q6) and z < 0.74 dAhffl-'QOian et al.. 2006b; Albert et aU2007a). 



Table 2. Measured velocity dispersions and resulting estimated black hole masses for the VHE blazars. 



Object 



o" [km s ' ] 
Reference 1 



log(M./Mo) 



cr[kms '] 
Reference 2 



log(M./Mo) 



o" [km s ' ] 
Reference 3 



log(M./Mo) 



log(M./Mo) 
averaged 



Mkn 421 


219±11 


8.29 ±0.18 


236 ±10 


8.42 


±0.15 


324 ±18 


8.97 ±0.12 


8.56 


Mkn 501 


372 ±18 


9.21 ±0.11 


291 ± 13 


8.78 ±0.11 






9.00 


lES 2344+514 


294 ± 24 


8.80 ±0.15 






389 ± 20 


9.29 ±0.12 


9.04 


Mkn 180 


209 ± 1 1 


8.20 ±0.19 


244 ± 10 


8.47 ±0.14 


251 ± 16 


8.52 ±0.15 


8.40 


lES 1959+650 






195 ± 15 


8.08 ± 0.23 


219 ±15 


8.28 ±0.19 


8.18 


BL Lacertae 












245 ± 16 


8.48 ±0.15 


8.48 


PKS 0548-322 


202 ± 24 


8.14 ±0.24 












8.14 


PKS 2005-489 












257 ± 16 


8.57 ±0.14 




H 1426+428 












269 ± 16 


8.65 ±0.13 




H 2356-309 












195 ± 14 


8.08 ± 0.23 




lES 0229+200 












363 ± 19 


9.16±0.11 




lES 1218+304 












191 ± 14 


8.04 ± 0.24 




lES 0347-121 












214 ±15 


8.24 ±0.19 




lES 1011+496 












219 ±15 


8.28 ±0.19 





References: (l) lBarth. Ho & Sargenll i2003l) : (2) lFalomo. Kotilainen & Trevej i2002h ; (3)IWu. Liu & ZhaneH2O03) . The velocity 
dispe rsions cr w ere translated into estimated BH masses using the M, — cr relation from lTremaine et al.l 120021). The cr values taken 
from lWu et al.l J2002) were indirectly determined using the fundamental plane of radio galaxies JBettoni et aU '2001). BH masses 
given in units of the solar mass, Mq . If more than one M, value is given, the average M. is used. Due to the possible different 
systematic errors of the individual data sets, the largest error was assumed as error of the average M, . 



using the virial BH mass estimate of IVIcLure & Dunlog ( l2002h and 
is given as log(M./]VIo) = 8.912 by Gu et al. (20oil) . For two of 
the blazars under study, PG 1553+113 and lES 1101-232, no M. 
estimations exist yet. 

Recent estimations of the BH masses for 452 AGN find them 
distributed over a large range of (10* -7 x 10') Mq wit h no evidence 
for depe ndencies on the radio loudness of the objects JWoo & Urrvl 
l2002al |b h. A recent study of the BH mass distribution of 66 BL Lac 
obiects lwoo et alj200a) reports an M. range of ( 10' -4 x 10') Mq 
and could also not find a correlation of M, with radio or X-ray 
luminosity. (In Sect. 15.51 the distribution of VHE blazars in lu- 
minosity and M. is discussed). As Fig. [3] shows, there is no de- 
pendence of the BH masses of the VHE blazars on their redshift, 
but they are rather flatly distributed in their BH masses between 
(10* - 10''^)Mo. Although AGN harbour BHs with M. > lO^^M©, 
up to now only blazars with rather massive BHs, M. > IO'^Mq, 
have been discovered in VHE y-rays, raising the question whether 



a physics reason is responsible for the non-detection of blazars with 
less massive BHs in the mass range (10' - 1O*)M0. There ex- 
ist studies that find radio-loud AGN, and there fore also b lazars, 
to be associated with BHs with M. > IO'Mq JLaoilbOOOl ) or at 
least on average with more ma ssive BHs than radio-quiet AGN 
jJVIetcalf & Magliocchettil 120061) . The latter authors also report a 
threshold BH mass for the onset of radio activity, with very little 
dependence of the radio output on the BH mass once above the 
threshold mass. Whether such a mass threshold is also at work 
for the VHE emission, remains subject for further studies at this 
point. The BH masses of the VFIE b l azars a re compared to those of 
375 AGN collected bv lWoo & Uffvl J2002ah in Fig.H The confine- 
ment of Seyfert galaxy measurements to low redshifts presumably 
is due to a selection effect: These are spiral galaxies and therefore 
expected to harbour comparatively low-mass BHs. Distant Seyfert 
galaxies (z > 1.0) might just not be luminous enough to obtain M. 
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Figure 3. Redshift vs. M, distribution for the known VHE y-ray emitting AGN. Tiie superimposed liistogram shows the BH mass distribution of the 375 AGN 
given in Fig.|4](linear vertical scale independent of z). 
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measurements. Conversely, quasars are too rare as to be found in 
small volumes and thus at small distances. 



4 INTRINSIC VHE y-RAY EMISSION PARAMETERS 

VHE y-ray observations enable us to look deep into the emission 
regions of blazar jets and thus convey information on the respon- 
sible particle acceleration and cooling processes. Here I study pri- 
marily the differential energy spectra in the VHE domain, which 
are summarised in Tab. [3] For the sources Mkn 421, Mkn 501, 
PKS 2155-304, lES 1959-1-650, and lES 2344-1-514 observations of 
clearly distinct flux states exist. Accordingly, for each of those 
sources two spectra, one 'low-state' and one 'high-state' spectrum, 
are considered. Low-state spectra are characterised by the absence 
of high (> 0.2 Crab unito flux levels and short-term variabil- 
ity (probably beyond instrumental sensitivity though), while flare 
spectra were obtained during outbursts of the respective sources 
(viz. the 1995 December 20 flare of lES2344-t-514, the 1997 flare 



' The Crab nebula exhibits a strong, constant VHE y-ray flux and is there- 
fore often considered a standard candle in VHE y-ray astronomy 



of Mkn 501, the 2002 flare of lES 1959-t-650 and the 2006 July 28 
flare of PKS 2155-304). At present for none of these sources, even 
for those with only one flux state, can a true baseline flux state be 
claimed, although low-flux states have been observed f Alb ert et al.l 
|2006a, 2007b). Long-term m onitoring cam paign s are currently per- 
formed to ad dress this issue JGoebel et al.ll2007l ; ISteele et alj2007l : 
|Puncdl2007h . 

The measured spectra suffer yy absorption on photons of the 
EBL as shown in S ect. [2l The i ntrinsic source spectra are recon- 
structed employing (lMazinl2003h the EBL 'low-lR' model given in 
iKneiske et alJ(l2004h . which assumes the least possible infrared star 
formation rate as allowed by galaxy counts and which is in reason- 
able agreement with other models (Prim ack, Bullock & Somervillj 
I2OO5I : IStecker. Malkan & Scullvll2006i) . Note that due to the fact 
that the VHE y-rays are attenuated exponentially with the optical 
depth, an accurate knowledge of the EBL is crucial for the individ- 
ual interpretation of the intrinsic VHE y-ray spectra. 

In the following, 1 will use two ob servables to chara cterise the 
VHE y-ray emission: The A'-corrected JHogg et alj2002l) luminos- 
ity at 500 GeV, VyLy = 47Tdl-(500GeWfF{500Gsy/(l+z))/(l+z) 
with the luminosity distance ^l and the intrinsic photon index 
r in the region around 500 GeV, which is determined by fitting 
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Table 3. Measured VHE blazar spectra, reconstructed intrinsic spectral indices and luminosities. 



Object 



Measured Energy Spectrum dF/dE 
[TeV^'cm-^s-'] 



Reference 



Intrinsic 
Slope r 



VyLy ^ 

[ergs"' sr"'] 



Mkn421 
Mkn501 
lES 2344+5 14 
MknlSO 
lES 1959+650 
BL Lacertae 
PKS 0548-322 
PKS 2005-489 
PKS 2155-304 



H 1426+428 

lES 0229+200 
H 2356-309 
lES 1218+304 
lES 1101-232 
lES 0347-121 
lES 1011+496 
PG 1553+1 13" 
PG 1553+113* 



(12.1 ± 0.5)10"'2(£/1.0TeVr3"9±'">'' 
(8.4 ± 0.5)10-i2(£/i.OTeV)-2'6+o.08 

(1.2±0.2)10-"(£/0.5TeVr2'5±'''2 
(4.5 ± 1.8)10-"(£/0.3TeVr3-3±''-7 
(3.4 ± 0.5)10-'2(£/1.0TeVr2''2±0'4 
(1.9±0.5)10-"(£/0.3TeVr3'i4±o.54 

(1.9 ± 0.4)10-'\£/1.0TeVr2-8*''-3 
(1.9 ± 0.7)10-'3(£/1.0TeVr'"'*''* 
(1.96 ± 0.12)10"'2(£/1.0TeV) 



-3.32+O.OS 



Q+0.46 
-0.27 



3.15 



+0.10 
-0.12 



(2.4+»-^)10-i2(0.7±0.2)<"'': 
x(£/1.0TeV)"^''''-o27 for E > 700 GeV 



for£<700GeV, 

) 



-2.6±0.6 



(2.9 ± l.l)10-"(£/0.43TeV) 

(2.34±0.37)10""'(£/3.0TeVr2 5±''2 
(3.08 ± 0.75)10~'^(£'/1.0TeV)-3''6±ot 



(8.1 ±2.1)10""(£'/0.25TeVr3"±"'' 
(4.44 ± 0.74)10-"(£/1.0TeV)-2**8±oi4 
(4.52 ± 0.85)10"'\£/1.0TeV)-3"'±°-23 

(2.0 ± 0.1)10-"'(£/0.2TeVr'"'±''5 
(1.8 ± 0.3)10- "'(£/0.2TeVr''2i±o.25 

(1.8 ± 0.3)10-"'(£/0.2TeVr''-'±''-25 



Aharonian et al. (1999b) 
Aharonian et al. (2001b) 
Albert et al. (2007b) 
Albert et al. (2006c) 
Albe rt et al. (2006a) 
Albert et al. (2007d ) 
Superina et al. (2007) 
Aharonian et al. (2005b'l 
^h^otuM^e^aL Q005a) 



Horan & Finlev (20011 . 
Aharonian et al. (2002) 
Raue et al. (20071 
Aharonian et al. (2006a) 
Albert et al. (2006b) 
Aharonian et al. (2006aJ 
Aharonian et al. (2007c) 
Albe rt et al. (IQOTe) 
Albert et al. (2007a1 
Albert et al. (2007a) 



2.85 + 0.58 
2.49 ± 0.84 
2.67 + 0.21 
3.06 ± 0.50 

2.37 ± 0.29 
3.17 + 0.25 

2.38 + 0.28 
3.52 ± 0.27 
2.43 ± 0.64 



1.39 + 0.30 
1.77 + 0.37 
1.97 + 0.40 

1.33 + 0.37 
1.76 + 0.14 
2.56 + 0.29 
3.68 + 0.68 

2.34 + 0.46 



(9.68 
(6.89 
(2.80 
(2.04 
(5.95 
(2.35 
(9.66 
(2.53 
(6.32 



+ 0.40) 
+ 0.41) 
+ 0.47) 
+ 0.81) 
+ 0.88) 
+ 0.62) 
+ 2.05) 
+ 0.93) 
+ 0.39) 



xlO-*' 
xlO« 
xlO« 
xlO« 
xlO'*3 
xlO« 
xlO' 



,42 



XlO' 
xlO' 



1.58 + 0.23 (7.42 + 2.81) XlO' 



,44 



(6.13 
(3.29 
(1.26 
(3.19 
(3.19 
(15.4 
(6.40 
(3.14 



+ 0.98) 
+ 0.80) 
+ 0.33) 
+ 0.53) 
+ 0.60) 
+ 0.77) 
+ 1.07) 
+ 0.52) 



XlO"' 
XlO"* 
xlO''5 
XlO'*'* 
XlO'*'* 
xlO"'* 
xW^ 
xl0'*5 



Mkn421'' 
Mkn501^ 
lES 2344+514^ 
lES 1959+650^ 
PKS 2155-304'' 



(23.40 + 0.73)10-"(£/l-0TeV) 



-2.32+0.03 

(2.50 + 0.16)10""'(£/1.0TeV)-222±o.04 

(5.1 + 1.0)10-"(£/1.0TeV)-2^''±°" 

(1.23 + 0.25)10""'(£'/1.0TeV)-2'8±0-l2 

(2.06 + 0.16)10"'°(£'/1.0TeV)-"i±o"6 for E < 340 GeV, 

(2.06 + 0.1 6) 10- '"(0.430 + o.022)(3«±oo5)-(2.7i+o.06) 

x(£:/1.0TeV)-3 "±005 for E > 340 GeV 



Krennrich et al. (20021 
Aharonian et al. (1999a) 
Schroedter et al. (20051 
Daniel et al. (2005) 
Aharonian et al. (2007b1 



2.09 + 0.30 
1.95 + 0.41 
2.20 + 0.31 
2.43 + 0.29 
2.28 + 0.40 



(1.08 + 0.031x10''' 
(1.39 + 0.09) X 10"' 
(6.65 + 1.30) XlO'*" 
(2.25 + 0.46) X 10"' 
(2.74 + 0.17) X 10"* 



r denotes the reconstructed (intrinsic) VHE spe ctral power-law index at 500 GeV and VyLy represents the source luminosity at 500 GeV. Both values were 



calculated from the measured spectra assuming a lKneiske et al.l 'low-IR" EBL density. " at an assumed z - 
during a flare state of the respective blazar. 



0. 1 . at an assumed z = 0. 3 . '^ spectrum measured 



the intrinsic spectra with pure power-laws of the form dF/dE = 
/o ■ (E/EoY^. These two parameters act as proxies for the peak 
position and the spectral shape on the falling edge of the high- 
energy bump, which cannot, as explained before, easily be deter- 
mined from the existing VHE data. For extraction of the luminosity 
and of the spectral slope, the region around 500 GeV was chosen 
because all blazars under study have measured spectra in this en- 
ergy region. All calculations and fits have been performed in en- 
ergy ranges where y-ray spectra for the respective blazars have ac- 
tually been observed, so that no extrapolations in energy regions 
not covered by the data were required. All in all, extragalactic 
source observations included in this paper cover the energy range 
85 GeV < £ < 1 1 Te V. For the determination o f the luminos i ty dis- 
tances the cosmological parameters given in Spere el et alj feOOTD 
were used: Cl^nh' 
ble constant Hq = 



= 0.127t°?^nb/^- = 
lOO-Zikms-'^"— -' 



0.0223!;;;™^ with the Hub- 



'Mpc 



:73+^kms- 



Mpc 



Tab.|3]also shows the resulting source luminosities in the VHE 
region and the spectral slope of the reconstructed source-intrinsic 



spectra. While the luminosities range from 



10"3 



ergs 
3 X 10"' ergs"' sr"' (w 10"' ergs"' sr"' to a; 3 x 10"* ergs"' sr 



sr 



for blazars in outburst), the photon indices of the reconstructed in- 
trinsic spectra vary between F = 1.4 - 3.3, except for lES 1101- 
232, which probably has an i ntrinsic spectrum peaking far beyond 
E = 1 TeV dAharonian et alj|2007a) . The VHE measurements lie 
close to, but generally above the maximum of the high-energy 



bump (which occurs at F = 2). Spectra with F < 1.5 are dif- 
ficult to explain in current ac celeratio n models |Malkov & F)niryl 
200 ll; lAharonian et al.ll2006J: but see iKatarzvnski et alj|2006dlbl; 
Stecker et al.ll2007 for models that would explain harder spectra). 
The rather hard intrinsic slopes found for some sources go in line 
with indications that the E BL absorption effects ar e still smaller 
than currently modelled (e.e. lAharonian et al.l2006al) : a lower EBL 
level would soften the intrinsic spectra inferred here, i.e. increase 
the value of F. 

Before instruments like MAGIC and H.E.S.S. became op- 
erational, the average observed photon index in BL Lac objects 
was F X 2.3. This raised the expectation that AGN would in 
general exhibit rather hard spectra in the VHE range, which also 
would be compatible with the average EGRET blazar spectrum 
(at MeV-GeV en ergies), found to have a slope of Feqret = 2.27 
jVenters & Pavli dou 2007). In fact, the average spectral slope at 
VHE has not changed much, although the scatter of F increased 
as a new population of objects with intrinsically rather soft spectra 
(Mkn 180, PKS 2005-489) has been tapped, and at the same time 
distant, hard-spectrum sources were found. 



For PG 1553+113 with its unknown redshift (see Sect. l5.6l for 
details), two possible distances (z = 0.1, z = 0.3) were assumed 
in this paper. The resulting 'intrinsic spectra', however, are only 
given for illustrative purposes and are not used for any conclusions 
throughout this study unless stated otherwise. 
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Figure 5. VHE y-ray luminosity VyLy vs. X-ray luminosity at 1 keV, vxLx for 17 VHE blazars. The two data points for PG 1553+1 13 (open crosses) are for 
assumed redshifts of z = 0.1 (v^Ly = 6.4 X W*^ ergsr"' s"') and z = 0.3 (VyLy = 3.1 X 10*^ ergsr"' s"'), respectively, and are not used in the fit and for 
determining the correlation coefficient (see text). It should be noted that VHE data points have an additional systematic error of typically 35 per cent. The 
systematic error of the X-ray luminosities is unfortunately unknown. 



5 CORRELATION STUDIES 

5.1 Correlation of X-ray, optical, radio and VHE y-ray 
luminosity 

In SSC models, the X-ray and the VHE emission are closely con- 
nected, owing to their common origin. While in some blazars clear 
evid ence for a corresponding correlation has been observed (Mkn 
421.lKrawczvnski et alj200ll : lBlazeiowski et al.l2005l : lAlbert et all 



2007d). the conne ction is only weak for other ones (Mkn 501, 
Albert et al .11 200711) or even non-existi ng during some flar e states 
(the lES 1959-1-650 orphan flare case. (Paniel et al.ll2005h . Fig.|5] 
shows VyL^ versus the X-ray luminosity at 1 keV (vxLx\ from 
ICostamante & Ghise llini 2002). Note that high thermal contribu- 
tions at 1 keV are unlikely and would imply a very high amount 
of gas and pressure. Unfortunately, the X-ray and VHE data have 
not been taken simultaneously. While VHE measurements during 
outbursts were not used in the Figure, variations in the X-ray do- 
main for the blaz ars under study are, a ccording to the compilation 
of X-ray fluxes in lDonato et al.l J200lh . not larger than a factor 4.6 
(most extreme object: Mkn 501) with an average of a factor 1.4 
and a variance of 2.2. According to expectations, a trend towards 
a correlation is visible. When including all data points except for 
those representing PG 1553-1-113, 1 find a correlation coefficient of 
r = 0.76^Q j', which is within 3.6 standard deviations different from 
zero. A linear fit to the data yields a slope of ;n = 1.11 ± 0.09 
OrL = 61.8/14). 

Recently, optical trigg ers lead to the succ essful discoveries of 
the VHE blazars Mk n 180 JAlbert et alj2006cl) and lES lOllH-496 
d Albert et alJl2007en . At times of lower optical emission the lat- 
ter in fact showed a lower fl ux and thus only yielded a marginal 
detection d Albert et al.l 120 08). This seems to imply the existence 
of a VHE - optical connection, which would be a convenient 
proxy for finding new VHE blazars. However, in the past no 



JBlazeiowski et alj|2005l: iRebillot et al.ll2006l ; JATbert et al.ll20073) 
or only weak dSucklev et alJ 1 19960 evidence for optical or ra- 
dio to VHE correlations was found for individual observations of 
Mkn 421 and PG 1553-1-1 13. 

The optical flux data at 5 500 A used here were collected f rom 
several source catalogues by ICostamante & Ghisellinil (2003), as 
were the radio flux data at 5 GHz. I converted them into luminosi- 
ties taking into account the appropriate luminosity distances d]^. 
Fig.|6]shows the corresponding correlations: the data in the VHE - 
optical plane feature a larger scatter than that in the VHE - X-ray 
data, while in the VHE - radio plane no clear trend is seen. 

VHE blazars might populate only restricted ranges in X-ray, 
optical, or radio luminosity distributions. To test this, I use X-ray, 
optical and radio data from the fu ll set of 246 sources considered by 
ICostamante & Ghisellinil d2002i) . I rejected the sources for which 
no redshift is known and converted the remaining 183 fluxes into 
luminosities. These blazars are compared to the VHE blazars in 
Fig. [7] In none of the three distributions can substantial deviations 
of the VHE blazars from the overall set of blazars be found. 



5.2 Correlations of intrinsic photon index and VHE 
luminosity with the synchrotron peak frequency 

In SSC models the VHE peak, identified with the inverse Comp- 
ton (IC) peak with a maximum at v'^^j^, resembles the form (e.g. 

iFossati et al.l 119981) of the synchrotron peak at v^^, displaced 



by th e squared Lorentz factor v|5..i,/ 



,sy 



-I 



dTavecchio et all 



1 19981) . iNie ppola. Tomikoski & Valtaoial ( l2006l) collected (non- 
simultaneous) multiwavelength data for a large (> 300) set of 
blazars, which includes eleven of the VHE blazars under study 
here. The data, covering frequenci es from the 5 GH z radio to the 2.4 
keV X-ray domain, were used bv INieppola et al.l to reconstruct the 
SEDs of these blazars and to infer the synchrotron peak frequency 
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Figure 6. (a) Correlation of VHE luminosity VyLy and optical luminosity vqLo- (b) Correlation of VHE luminosity VyLy and radio luminosity vrLr. The two 



data points for PG 1553+113 (open crosses) are for assumed redshifts of z = 
respectively. 
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Figure 7. Histograms of (a) 1 keV X-ray, (b) 5500 A optical and (c) 5 GHz radio luminosity of the VHE blazars and all those blazars considered by 
ICostamante & Ghisellinil 420021) with determined redshifts ( 1 83 out of the total set of 246 blazars). 



Vpcak- I test for correlations between v ^^j^ and the VHE luminosity 
and spectral slope (Fig. [8}. According to expectations from SSC 
models, I find a correlation of the photon index F with v J^. The 
correlation coefficient is r = -0.63^Q^g and a linear fit to the data 
as shown in Fig.[8ji yields ax^^^ = 12.1/8. Fig.[8j> shows the corre- 
sponding data points in the VyLy - v ^^j^ plane, in which no correla- 
tion is apparent. Finally, in Fig.|9l I check whether the LBL-HBL 
transition of the VHE blazars is connected with the BH masses of 
their host galaxies. Ten of the VHE blazars, for which both v ^^j^ 
and M. measurements are available, show no trend towards a cor- 
relation. 



5.3 Correlation between the intrinsic photon index and the 
VHE y-ray luminosity 

Fig.llOti relates the intrinsic photon indices F to the VHE y-ray lu- 
minosities VyLy. For all objects under study (excluding flare states 
and the PG 1553+113 data), the F - VyLy correlation reads as 
F = Fo + m-log,o(v^L^)withFo = 38.5 ±6.8 and m = -0.82 + 0.15. 
The fit which led to this correlation has a x^^^^ = 40.18/14. The 
line at F = 2 denotes the spectral slope at which a maximum in 
the SED occurs. Looking at the spread of the observations in F as 
a function of VyL^, I notice that the distribution sharpens towards 
F = 2. This might reflect that the highest luminosity occurs at 
F = 2. Thus the spread of the data reflects the spread of the shapes 
of the high energy peaks. The general behaviour above F a; 2 - 
the higher the VHE y-ray luminosity, the harder the spectrum - 
can within SSC models be described with a moving IC peak to- 
wards higher energies with increasing luminosity. The five blazars 



with observed spectra at quiescent and flare states partially be- 
have i n a similar matmer (Fig. 1 10b): Mkn 421 I Krawczvnski et al.l 
l200ll : lB}azeiowski et alj|2005,). Mk n 501 ( lAlbert et alj|2007 J ) and 
lES 2344+514 jAlbert et alJl2007bh also show a spectral harden- 
ing during high-flux states, flares and outbursts. Fig. Illb shows 
the corresponding luminosity differences A(VyLy) and slope differ- 
ences AF. Mkn 501 and lES 2344+514 show a similar change in 
spectral slope and a dynamical range of A(VyLy) x 20. The lu- 
minosity increase of Mkn 421 observed up to now is much lower 
with A{VyLy) » 10. While the spectral slope of Mkn 421 also 
hardens with increasing luminosity, lES 1959+650 and PKS 2155- 
304 show within errors no variation in their spectral slopes during 
flares, while their luminosities increase rather drastically by a fac- 
tor of 40 and 50, respectively. When plotting the luminosity differ- 
ence versus the BH masses (Fig.lllb). a trend, although broken by 
lES 1959+650, towards blazars with more massive BHs showing 
higher dynamical ranges of their emission, is found. 



5.4 A test for selection efTects: redshift dependencies? 



pected, but may identify selection effects in the data set and/or 
an inaccurate EBL model. Conspicuously, very hard (F <K 2.0) 
spectra have up to now been only reconstructed for rather dis- 
tant (z > 0.1) blazars (Fig. I12h) , a trend also observed by 
IStecker, Baring & Summerlinl ( 120071) , who used a smaller set of 
VHE blazars. At the same time, none of the measured spectra of 
nearby sources shows F much smaller than 2.0, although for these 
blazars no strong EBL modifications apply and the measured spec- 
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Figure 8. (a) Photon index F vs. synchrotron peak frequency v , and (b) VHE kimin osity vs. synch r otron peak frequency v , for eleven VHE blazars. 



The synchrotron peak frequencies and the LBL-IBL-HBL classification are taken from Nieppola et alj 420061) . The two data points for PG 1553+113 (open 
crosses) are for assumed redshifts of z = 0.1 (E = 3.68 ± 0.68; v^Ly = 6.4 X lO''^ ergsr"' s"') and z = 0.3 (E = 2.34 ± 0.46; VyLy = 3.1 X lO^^ergsr"' s"'), 
respectively, and are not used in the fit and for determining the correlation coefficient in the left Eigure (a). 
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Figure 9. Estimated BH mass vs. synchrotron peak frequency v ^ . for those VHE blazars with both quantities known. 



tra should not differ substantially from the intrinsic ones in the en- 
ergy region studied here. Why are only blazars with rather hard 
intrinsic spectra visible at large distances (z > 0.1)? Soft spectra 
certainly fall more easily below the current instrumental sensitivity 
limits. Another explanation for the prevalent hard spectra at large z 
is an overcorrection of the EBL attenuation effects. 

Fig. 112b shows the corresponding distribution of VHE lumi- 
nosities VyLy as a function o f z- Two curves indicate the sensitivity 
limit of current lACTs (e.g. lAharonian et al.ll2006br) . a significant 
detection of 1 per cent of the flux of a Crab nebula-like source in 
25 hours, and the sensitivity of previous (before 2002) lACTs of 
w 10 per cent of this flux. Interestingly, most of the sources found 
at z > 0.05 seem to be rather low-luminosity sources in the sense 
that their luminosity is not much higher than the current instrumen- 
tal sensitivity allows for. This means that not only the substantially 
lower energy thresholds of the current lACTs (< 100 GeV as com- 
pared to X 300 GeV until 2004), but also their increased sensitivity 
enabled some of the new blazar discoveries. Exceptions to these 
detections close to the current sensitivity limit are (trivially) the six 
VHE blazars discovered before 2002, lES 1011-1-496 (but discov- 
ered during an optical flare, while it s low-flux state seems to l ie be- 
low the current instrumental reach. I Albert et alj|2007el . 120080 . and 
lES 1218-1-304, which in fact seems to have been detected in a non- 
flare state. Not e further th at H 1426-1-428 has not yet been detected 
after 2002 (e.g. lKrawczvnski 2007 : Albert et al. 2008) . which also 
currently places it below sensitivity limits. 



5.5 Correlations of intrinsic VHE y-ray emission parameters 
with the black hole mass 

The y-ray production is thought to take place at shock fronts inside 
the AGN jets at very close (sub-parsec ) distances from the central 
BH JJester et al.l2006l : IUchivama et al.l 2006). While the jet produc- 
tion and collimation mechanism is still elusive, accepted models 
are generally based on m agnetohydrodynamic (Blandford & Paviig 
ll982l : lKudoh et alJI 19991) or electromagnetic jet models. In the lat- 
ter, a Poynting flux domin ated flow is launched from a Kerr BH 
(IBlandford & Znaiekll 19770 or from the accretion disc I IBlandfordl 
1976). The conversion from Poynting dominance into particle dom- 
inance is not yet understood. The properties of the blazar y-ray 
emission are expected to be connected to the properties of the cen- 
tral BH, like M. and the BH spin, since scaling laws govern BH 
physics (McHardv et al. 2006), in part i cular length an d time-scales 
of flows ( Mirabel & Rodriguez! iT999l : iMirabell 120041) . e.g. the or- 
bital period of the last stable BH orbit. Currently, only M. can 
be reliably estimated; the BH spin remains inaccessible by large. 
Moreover, the environment in which the BH is embedded might be 
equally important; one of its properties, the ac cretion rat e, is indi- 
rectly accessible through the (radi o) jet power JLiu et al. 2006), or 
from multi wavelength modelling JMaraschi & Tavecchia.2003.) . A 
previous study of the connection of spectra l properties and M. for 
five VHE blazars JKrawczynski et al.ll2004l) did not find any corre- 
lations with the BH mass. 

For 15 VHE blazars with known BH masses, I neither find 
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a correlation between M. and the spectral slope F (Fig. 113b). nor 
between M, and the VHE y-ray luminosity (Fig. 1 13b). While the 
VHE blazar set disfavours a dependency of the VHE y-ray emission 
properties studied here on M., the uncertainties of the M. determi- 
nation are still rather large and might conceal otherwise interesting 
physics. In particular, indirectly inferred a values from the funda- 
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Figure 11. (a) Evolution of intrinsic spectral index F and source lumi- 
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mental plane or bulge luminosity measurements have rather large 
uncertainties. To improve on the M, uncertainties, it would be de- 
sirable to obtain direct measurements of a for all VHE blazars. 
Moreover, VHE emission properties may also depend sensitively 
on the BH spin, the accretion rate, or on properties of the accel- 
eration region in the jet. Also results on timing properties (see 
Sect. l5.7b support such claims. 



5.6 An upper distance limit for PG 1553-hll3 

The redshift determination for blazars is challenging, as these AGN 
generally exhibit only weak spectral lines. In several attempts, no 
emission or absorption lines could be found in the optical/TR spec- 
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trum of PG 1553+113; see e.g. IWagner et aO ( l2007h for results of 
a recent Very Large Telescope (VLT) campaign. T he frequently 
cited initial determination of z = 0.36 JMiller & Gre en 1983) was 
found to be based on a misidentified emission line and could not be 
repro duced JFalomo & Treveslll99d ; Falomo, Scarpa & Bersanellil 
Il994) . VLT optical spectroscopy JSbarufatti et al.ll2006l) yields a 
lower limit of z > 0.09, while the analysis of Hubble Space Tele- 
scope im ages leads to the predic tion of a redshift in the range of z = 
0.3 - 0.4 I ITreves. Falomo & Us lenghi 200 71). Indirect methods em- 
ployin g the maximum slope JAharonian et alj|2006bl ; lAlbert et al.l 
l2007ah or the shape JMazin & Goebell |2007|) of the VHE y-ray 
spectrum find upper limits of z < 0.74 and z < 0.42, respectively. 

With increasing distance, the luminosity of PG 1553+113 has 
to increase stronger than quadratic due to EBL yy absorption as 
to sustain the measured VHE flux. Fig. [14] shows the source lu- 
minosity (1) when only considering the distance effect and (2) 
when also taking into account the EBL absorption. Due to the 
exponential behaviour of the EBL attenuation, the latter effect is 
by far dominant. I assume here that PG 1553+113 is an 'off the 
shelf blazar, i.e. with no extraordinarily high luminosity VyLy. 
This assumption is difficult to quantify; given the overall dynam- 
ical range of the (non-flare) blazar luminosities in this study of 
w 75, 1 consider the case in which the luminosity of PG 1553+113 
is not more than 30 (1000) times higher than the highest luminos- 
ity found in the sample. Further, PG 1553+113 has in two years 
of observations not shown any apparent flaring b ehaviour; within 
a factor of three, the measured flux was constan t JAharonian et al.l 
l2006bl ; [Albert et al.ll2007l iBenbow et alj|2007h . Therefore I con- 
sider only sources in non-flare states, of which lES 1218+304 with 
VyLy = 1.3 X lO'*' erg s"' sr"' is the most luminous one. A 30-times 
higher luminosity then implies a limit of z < 0.45, while an ex- 
treme luminosity of VyLy = 1.3 x 10''^ergsr~' s"' yields a limit 
of z < 0.64. These limits do not only depend on a good knowl- 
edge of the EBL over a wide range in redshift, but also on the 
assumed maximum VHE blazar luminosity that strongly depends 
on the Doppler factor S. In any case, either a strikingly high lu- 
minosity or a very high S is needed to explain the observations 
should PG 1553+113 be more distant than z > 0.35. Note that 
^ < 20 suffices for most of the blazars modelled up to now dur- 
ing non-outburst time s, as also for flare observations (e.g. Mkn 421, 
iMaraschi et alj 19991). SSC modelling for PG 1553+1 13 res ulted in 
6 = 21 iCostamante & Ghisellinill2002l ; [Albert etalj2007ah . 

The difficulties in finding emission and absorption lines might 
indicate a very close alignment of the jet axis of PG 1553+113 
to our line of sight. Very-long baseline interferometry (VLBI) 



imaging is available for some VHE HBLs detected before 2002 
jPiner et al.[[T99^ ; [Edwards & Piner"2002'; 'Piner & Edwards[ [2004h 
and for BL Lacertae (JDenn, Mutel & Marscher 2000), while addi- 
tional jTieasuiCTients_oftherecently found VHE blazars are under- 
way dPiner. Pant & Edwards[[2007,) . VLBI essentially confirms the 
close aligmnent of the jet to our line of sight on sub-parsec scale 
and finds opening angles of typically few degrees. These results 
cannot yet be used for quantitative correlation studies though. 

We have to eventually consider the possibilities that 
PG 1553+1 13 is a rather distant source and the blazar populations 
at large distances show significantly different properties than the 
close-by objects at z < 0.2, and that such very extreme objects are 
so rare that a sufficiently large volume had to be probed to find one 
of them. 



5.7 Correlations with the X-ray duty cycle and the VHE 
variability time-scale 

Following a method described in [Krawczvnski et al.[ yOOJ), I de- 
termine the time fraction ('duty cycle') for which the (2 - 10) keV 
X-ray flux exceeds the average flux by 50 per cent. In this paper, 
for a blazar to be regarded 'on duty' this deviation is additionally 
required to be significant on the 3o" level. 2-10 keV X-ray light 
curves are obtained from the All-Sky Monitor detectoQ on board 
the Rossi X-ray Timing Explorer (RXTE) and are available from 
1996 January 5 on. Fig. [15] shows the corresponding light curves 
along with the resulting duty cycles. Objects which are classified as 
extreme BL Lacs (Mkn 501, lES 2344+514, H 1426+428, H 2356- 
309 and PKS 0548-322. lc"ostamante et alj200l[) show substantially 
higher activity than the other blazars. Note also the outstanding role 
of Mkn 421 with its substantially higher duty cycle as compared to 
all other objects. Fig. 116b shows the duty cycle as a function of 
the BH mass for the VHE blazars. No trends are visible, except 
for the observation that the three blazars with the most massive 
BHs (log(M./MQ) > 8.8) do not have duty cycles in excess of 17 
per cent. This goes in line with a speculation of an anticorrelation 
between the X-r ay flare duty cycle and M . seen in the five VHE 
blazars studied in [ Krawczvnski et al.[ ( [2004[) . The distribution of the 
X-ray duty cycle as a function of luminosity is found to be rather 
flat (Fig. 116b). supporting the claim that variability on all scales is 
a defining property of blazars. Note, however, that Mkn 501 and 
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PKS 2155-304, the objects which show the fastest variability in 
VHE y-rays, are among the objects with a rather low duty cycle. 

Blazars are characterised by a highly variable emission. In 
particular the VHE y-ray emission is often found to be more vari- 
able than the emission at other wavelenghts. Still, comprehensive 
data on the temporal behaviour of VHE blazar emission has of- 
ten been collected only for blazars in outburst, and thus, over 
short time spans. Although VHE light curves ranging over ten 
years' wo rth of measurements were collected occasionally (e.g. 
[Albert et al. 2007b , f; Tluczykont et al. 2007) , the sampling is only 
sparse, and continuous, unbiased l ong-term monitoring campaigns 
have not been star ted until 2005 JGoebel et al.ll2007l : 1steele et all 
l2007l : |Punchll2007l) . 

Variability time-scales t are intimately linked to the exten- 
sion of the region R from which the observe d em i ssion originates 
by the causality condition R < (5cr/(l -I- 7). ICuil d2004) has sug- 
gested that the flare hierarchy seen in Mkn 501 long-term data 
implies a scale-invariant nature of the flare process and that there 
might not be any fundamental difference among long, intermediate 
and rapid flares. Although not thoroughly understood, the flares in 
blazars might be related to internal shocks in the jet CRee s 1978a; 
ISpada et al.ll200lh or to major eje ction events of new components 
of relativistic plasma into the jet JBottcher. Mause & Schlickeiseil 
I1997I : iMastichiadis & Kirkiri997h . Different flare time-scales thus 
may be caused by a hierarchy of inhomogeneities in the jet, ener- 
gised so as to produce flares. As a timing property of the VHE y-ray 
emission, the minimum flux doubling times for the VHE blazars 
were collected from literature and plotted versus M. (Fig.ll7b) and 
the VHE luminosity (Fig. 117b). Because of the very limited data 
base of VHE variability measurements, which also is biased to- 
wards high-flux states and outbursts, time-scales are mostly given 
as upper limits, which disables strong conclusions. 

The observed VHE flux doubling times do not scale with the 
BH mass (Fig. 117b). which may simply mean that (1) the flaring 
mechanism is working in a much smaller region than the BH ra- 
dius/least stable orbit and more importantly (2) the BH and its prop- 
erties as such do not influence the flaring process substantially, and 
the embedding environment of the BH and the jet environment play 
more dominant roles (e.g., the accretion power of the system). The 
extremely short doubling times of r < 5 min recently found for 



Mkn 5 01 JAlbert et al.l2007 ^ and PKS 2155-304 JAharonian et"aD 
l2007bf) clearly do not support such a scaling either; even more, they 
are vastly incompatible with the rotation period of particles on the 
last stable orbit, which e.g. for the BH in Mkn 501 is T = 8.4 d. 
Additionally, the size scale implied by r < 5 min requires large 
Doppler factors S in the y-ray production regions as to avoid self- 
absorption. In contrast to the expected scaling behaviour of the flow 
properties around BH with their masses, the three AGN that host 
rather massive BHs, PKS 2155-304, Mkn 501 and Mkn 421 were 
found to exhibit the shortest variability time-scales. Until minute- 
scale flaring was found in PKS 2155-304, one could have argued 
that short time variability could only be measured for high fluxes 
due to the proximity of the respective sources: Mkn 421 and Mkn 
501 are the closest blazars at z < 0.034, while PKS 2155-304 is 
located at z = O.I 16. The present observations could, however, still 
be the cause of a selection effect: Minute-scale flares were found 
only during exceptional high-flux states of the respective sources 
so far. 

In Fig. [TS] I translate the BH masses into the corresponding 
gravitational radii r^ = GM/c^ (=0.5 Schwarzschild radius) and 
the observed minimal variability time-scales into the sizes of the 
emission regions R < 5ct/(I -i- z), assuming S = 10. For blazars 
in which fast variability has been observed, the extreme compact- 
ness of the emission region is apparent, being clearly comparable 
to or smaller than the Schwarzschild radius of the central BH. Even 
if the flares were driven by extremely large Doppler factors, say 
6 ~ 100, this would still result in emission regions smaller than the 
Schwarzschild radius for the Mkn 501 and PKS 2155-304 flares. 



6 CONCLUDING REMARKS 

Before the new generation of lACT became operational, only six 
firmly detected extragalactic VHE y-ray sources were known (e.g., 
Mori 2003); not for all of them differential energy spectra had 
been inferred. To date, there are 17 BL Lac objects and the FSRQ 
3C 279 known to emit VHE y-rays. This substantially enlarged set 
of blazars called for a synoptic study. I collected and derived in- 
trinsic properties of the VHE y-ray emission (luminosity, spectral 
hardness, temporal properties) and further included X-ray, optical 
and radio emission properties. As an accessible property of the BHs 
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Figure 15. X-ray (2-10 keV) light curves of VHE y-ray emitting blazars. The flare duty cycle, i.e. the fraction of time in which the respective object significantly 
exceeds its average X-ray flux by 50 per cent is shown. For convenience the ranges of the vertical axes are fixed for all plots. FSRQ: Flat spectrum radio quasar, 
LBL: low-frequency peaked BL Lac object; HBL: high-frequency peaked BL Lac object, synchrotron peak in the UV/X-ray range; XBL: extreme BL Lac 
object. 
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in the centres of the blazar host galaxies, the BH mass was also 
used. The studies yield the following results. 

(i) So far, only blazars with M, > 10^'Mq show VHE 7 emission. 
Whether this experimental finding constitutes a selection effect or 
will reveal interesting physics, remains to be seen. 

(ii) The VHE luminosity of the blazars under study and the corre- 
sponding X-ray luminosity show a hint of a correlation, as expected 
from leptonic acceleration/SSC models. 

(iii) A correlation between the spectral slope in the VHE region 
and the peak location of the synchrotron peak is found. Such a cor- 
relation is also expected from SSC models. 

(iv) There were no correlations found between the 7 emission 
properties and the M. of the galaxies that host the blazars. Also, 
no correlation of M. could be observed with the flare duty cycles 
and the flare time-scales. Thus, VHE 7-ray emission properties may 
not dominantly depend on M.. Other possibly interesting BH pa- 
rameters are not yet within instrumental reach. 

(v) There is an indication that the VHE 7 luminosity is correlated 
with the spectral hardness. This correlation can be formulated as a 
decrease of AF x 0.82 per decade of luminosity. 
(vi) This behaviour is also found for some individual blazars that 
were observed in different emission states: There are indications 
that in some variable sources the observed spectra become harder 
with increasing luminosity, while in others no hardening is found. 

Investigations of the temporal properties of the X-ray emission 
show that the blazars with the most massive BHs have rather low 
duty cycles. The temporal behaviour of the VHE 7-ray emission is 
not well studied for most of the VHE blazars, but the present data 
do not seem to support a scaling of the flux doubling time-scales 
with the BH masses. 

Obviously we are still dealing with a low number of sources 
and certainly with an incomplete source sample, which leaves re- 
gions in the parameter space empty. Nevertheless, the rather im- 
pressive number of 17 VHE blazars has permitted first comparative 
studies. It also allows conclusions concerning EBL models: A hint 
at a marginal correlation between the intrinsic spectral hardness and 
the source distance is likely due to a n EBL over-prediction. T his re- 
sult was also recently quantified bv lAharonian et alj ( l2006ar) using 
the two distant blazars lES 1101-232 and H 2356-309. In clarify- 
ing the situation, more distant sources particularly with intrinsi- 
cally hard spectra (as e.g. lES 0229+200) will play a decisive role; 
PG 1553+113 with its yet undetermined distance might also turn 
out to be a good candidate for such analyses once a measurement 
of its redshift succeeds. In the meantime, the observed luminosity 
distribution of the studied VHE 7 blazars was used to constrain the 
unknown distance of PG 1553+113 by assuming that the properties 
of this blazar are not too different from the most extreme objects 
in the blazar sample. Conversely, a large distance of PG 1553+113 
implies an unusually high luminosity or an unusually high Doppler 
factor. 

It can be said that the era of blazar astronomy has been entered 
- astronomy being understood as the study of generic properties of 
a given class of objects. With the currently known VHE blazars and 
the given instrumental sensitivity, VHE blazar astronomy starts to 
become less biased: Not only blazars with hard spectra or during 
outbursts have been detected, but also low (quiescent?) emission 
states in the VHE range are now being observed and studied. 

With continuing discoveries of new sources, it will become 
easier to disentangle EBL absorption effects and intrinsic absorp- 
tion effects from the measured spectra. Already currently, the VHE 
blazar sample contains groups of objects at very similar distances. 



e.g. lES 1218+304 and lES 1101-232 (Az = 0.004), or the triplet 
lES 2344+514, Mkn 180 and lES 1959+650 (Az = 0.003). Such 
groups are ideal for direct comparisons between the respective in- 
dividual objects, as they are subject to very similar EBL attenua- 
tion. 

Since 2003, on average three new blazars have been detected 
per year. When comparing this rate to the ten years it took to dis- 
cover the first six VHE blazars, one can with good reason hope to 
accumulate many more sources in the near future, which will refine 
synoptic blazar studies as the pioneering one conducted here. In 
addition, the first observation of new object classes like LBLs and 
FSRQs will certainly help to sharpen the VHE view of a possible 
blazar sequence and the underlying physics in the near future. 
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